Calculation of collective modes for the Bose-Hubbard model with confinement 
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The collective excitations in the Bose-Hubbard model in a trap are studied by means of numerical 
diagonalization in one dimension. The strength function is calculated for monopole and dipole 
perturbations, and moments of the strength function are utilized in order to obtain information 
about the collective behavior under external forces. In the superfluid regime, the spectrum is found 
to be exhausted by one single frequency, while in systems that contain a Mott insulating plateau 
several frequencies are excited. An explanation of recent experimental findings in terms of a Mott 
plateau is suggested. 



Recent years have seen a rapid development in the 
study of atoms confined in optical lattices. One of the 
virtues of such systems is that they provide an almost 
perfect realization of the bosonic Hubbard model. Cor- 
rections such as nearest-neighbor interactions are under 
control and can be made very small, and furthermore, 
the parameters of the model are widely tunable exper- 
imentally. In the case of spinless bosonic atoms the 
Bose-Hubbard model is known to support a Mott tran- 
sition, and it has also been demonstrated experimentally 
by Greiner et al. The addition of spatial confine- 
ment adds complexity to the well-known Bose-Hubbard 
model. In these inhomogeneous systems superfluid and 
Mott insulating phases may coexist in spatially separated 
regions 0- However, a good experimental diagnostics is 
still lacking for the detection of this interesting type of 
state. 

This paper wishes to forward the idea of using oscil- 
latory modes in order to discriminate between different 
types of many-body state for these systems. Exciting os- 
cillatory modes is an experimentally well-developed tech- 
nique for studying the many-body physics of trapped 
condensates in the absence of optical lattices and has 
recently been applied also to the case with an optical 
lattice present 0, _5] . Physical limits and basic depen- 
dences of the many-body eigenstates and energies were 
studied in an earlier paper |g (cf. Q)- Here we wish to 
relate the eigenvalue spectrum to the collective oscilla- 
tions of the system and pose the physical question, what 
happens when the system is perturbed? Just as in Ref. 
, a five-particle system shall be studied using exact di- 
agonalization, but the methods and concepts introduced 
here can be applied in numerical calculations using meth- 
ods such as Monte Carlo, density-matrix renormalization 
group or mean- field calculations. 

The system under consideration is a gas of bosonic 
atoms, confined in a harmonic trap superposed with an 
optical lattice formed by standing laser waves. By mak- 
ing the trap sufhciently tight in two directions, the sys- 
tem is effectively made one-dimensional, and when the 
optical-lattice potential is sufficiently strong the system 
is well described by the Bose-Hubbard Hamiltonian 



The units are here chosen so that the tunneling strength 
is unity. U is the on-site interaction strength, tOt is the 
trap frequency, and r is a dimensionless site index; the 
operators ar,al destroy and create a particle at the site r, 
respectively, and obey Bose commutation relations. The 
sum subscripted < rr' > runs over all pairs of nearest- 
neighbor sites as usual. 

We wish to know how the system responds to a per- 
turbation of the type 



F = y^v{r)alar, 



(2) 



where v{r) is a perturbing potential applied to the 
bosons. The strength function associated with the op- 
erator F is H 



5M=5^|(fc|F|0)p<5(. 



Ek)i 



(3) 



where k labels the many-body eigenstates of the system 
and Ek are the eigenenergies of the many-body system 
(in units where h = 1, and defining the ground-state 
energy as the zero level). 

In the present exact diagonalization study, the strength 
function can be obtained exactly. However, considering 
the outlook of using numerical methods where larger sys- 
tems can be treated, but the whole eigenvalue spectrum 
cannot easily be obtained, approximate methods must 
be considered. Reference describes a moment method 
(cf. that is well suited for describing the exci- 

tation spectrum for confined bosons in optical lattices. 
The moments of the strength function. 
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can, namely, be written as ground-state expectation val- 
ues of commutators involving the Hamiltonian and the 
operator F. For the three moments that shall be used 
here, mo, mi and m2, one obtains 

mo = {0\F^\0)-{0\F\Of, 
mi = 1{0\[F,[H,F]]\0), 



(1) 



1712 - mF,H][H,F]\0). 



(5) 
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The moment rUp has the dimension of frequency raised to 
the pth power, times the dimension of F^, and therefore 
the quotients between different moments define frequen- 
cies: from the ones hsted above we define ujiq = mi/mo, 
W20 = ('7i2/»Tio)^/^, and a'21 = m^/mi. In the case where 
only one eigenstate is excited by F, these moment-related 
frequencies coincide and are equal to the frequency of 
that state, but more generally, they take on values be- 
tween the frequencies of those modes that are appreciably 
excited by F , and the spread between them is related to 
the spread of the strength function. Observe that very lit- 
tle information about the individual excitation frequen- 
cies can be obtained from the moment calculation, ex- 
cept for the the fact that as p — + — 00, the frequencies 
LDp^p-i approach the lowest-lying eigenstate with a sym- 
metry consistent with the perturbation F The three 
frequencies are, however, average excitation frequencies, 
in a sense that is made precise by considering their defi- 
nition: 

and correspondingly for the other two frequencies. Fur- 
thermore, by writing toq as 

mo = ^|(fc|i^|0)^ (7) 

it is clear that the zeroth moment is related to the proba- 
bility for any excitation at all to happen, and it therefore 
gives information about the amplitude of the oscillations 
excited by the perturbation. 

Let us now explore in the setting of an exactly diago- 
nalizable system the relation between the strength func- 
tion and the moment-related frequencies ujpq. The details 
of the diagonalization scheme were explained in Ref. 0. 
The system consists of = 5 particles in one dimension 
subject to the Hamiltonian, Eq. (Q. The perturbing po- 
tential is assumed to be on the form of Eq. Q with the 
two choices of potential function 

v{r) — r , dipole perturbation (tilting), 
v{r) — , monopole perturbation (compression). (8) 

The exact strength function can be found by calculating 
the strength 5*^0 = |(fc|F|0)p for each excited state k. 
The strength function for the finite system is a sum of 
delta functions, and can of course not be plotted, but it 
is visualized in Fig. ^ by plotting a stem of height Sko 
at the position Ek of each excited state. For Fig. ^ the 
parameters are chosen as oJt = 0.3, U = 0.1, describing 
a weakly interacting system in a shallow trap, and the 
perturbing potential v is chosen as the monopole pertur- 
bation. It is seen that only two almost-degenerate modes 
are excited with appreciable strength. These are the two 
lowest-lying monopole modes, |2) and |3), at E2 « 0.5524 
and £'3 K, 0.5776; the excitation of the latter is weaker 
by approximately a factor 10. The strength for the third- 
strongest excited state is weaker again by another factor 




FIG. 1: Visualization of the strength function for a monopole 
excitation of a weakly interacting system, with ujt ~ 0.3, U = 
0.1. The height of the stems represent the strength Sko for 
each excited state \k), and their horizontal position indicates 
the frequency relative to the ground state of the excited state 
in question, u = Ek- The crosses mark the magnitude of three 
moment-related frequencies, in increasing order, ojio, 1^20 and 
CJ21. (However, the positions of the three crosses coincide 
completely in this figure.) The inset shows the density profile 
as circles and the quantum fluctuations in the local density 
as asterisks. All quantities are plotted in the dimensionless 
units employed throughout the paper. 

20. The strengths for the eigenstates that have dipole 
symmetry are suppressed by about ten orders of magni- 
tude, which corresponds to the numerical accuracy of the 
algorithm. This system is an almost pure Bose-Einstein 
condensate, its density resembling a discretized Gaus- 
sian, as seen in the inset of Fig. ^ Each of the sys- 
tem's collective modes is almost entirely exhausted by 
one many-body eigenstate or, as in the monopole case, 
an almost-degenerate pair of them. 

Included in Fig. ^ as crosses, are the three moment- 
related frequencies copq, which can be calculated using 
only ground-state properties. It is seen that all three fre- 
quencies lie close together - the symbols coincide almost 
completely so that only one is visible in the plot - and lie 
close to the main excitation frequency E2, as expected. 
The magnitude of the zeroth moment is toq = 30.6515; 
this number makes more sense when compared to values 
of mo for other parameter choices. It can be roughly in- 
terpreted to mean that application of the potential v{r), 
which means quadrupling the trap potential in this case, 
yields a population transfer rate of more than order unity 
to excited states and the oscillations should be clearly 
visible. 

Figure 13 displays the results for the dipole excitation 
using the same parameters LUt = 0.3, U = 0.1. The pic- 
ture is similar, but even simpler than in the monopole 
case: only one mode is appreciably excited, which in 
physical terms means that a perturbation with dipole 
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FIG. 2: Same as Fig.^ but here the perturbation has dipole 
symmetry. 
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FIG. 3: Same as Fig. Q but here the trap strength is ojt = 1 
and the coupling is intermediate, U = 20. In the inset, the 
asterisks represent the quantum fluctuations in the density 
magnified by a factor 10. 



symmetry results in excitation of the lowest eigenstate of 
the many-body system, and thus oscillation at one well- 
defined frequency. The zeroth moment is tuq = 1.9928, 
again of order unity. 

Another example is shown in Fig. 13 The parameter 
choice uJt = 1, U = 20 corresponds to a mixed state 
with a rudimentary Mott insulating plateau in the center. 
At the present there is very little understanding of the 
oscillatory modes in this type of regime and little hope 
for an analytical approach. The results for the strength 
function is displayed in the main panel of Fig.|31 Clearly, 
the strength is now spread out over many eigenstates, 
although the largest strength is about a factor 10 larger 
than the next-largest. In experiment one can expect to 
observe a monopole oscillation with several frequencies, if 



FIG. 4: Same as Fig. but here the trap and the coupling 
are stronger, uJt = i and U = 100. In the inset, the asterisks 
represent the quantum fiuctuations in the density magnified 
by a factor 100. 



damping is weak enough and the precision is good enough 
to resolve the different frequencies. The results of the 
moment calculation are likewise included in Fig. |3| The 
three moment-related frequencies now differ appreciably, 
as expected from the spread of the excitation strengths. 
As noted above, the moment calculation does not give 
direct information about the individual frequencies that 
are excited but the ujpg lie between the frequencies of the 
most strongly excited modes. The zeroth moment for the 
Mott-plateau system is found to be mg — 1.9596. This 
indicates that again, excitation should result in visible 
oscillations. 

The final example is in the limit of strong coupling in 
a tight trap, U — 100 and ujt = 4. This system is Mott 
insulating Figure 01 depicts the strength function and 
the moment-related frequencies for this system. Just as 
in the Mott-plateau case, the strength is distributed over 
several modes, which shows up as a difference between 
the moment-related frequencies. Observe, however, the 
magnitude of the strength function: correspondingly, the 
zeroth moment is now very small, mo = 0.0094. The rea- 
son is that the ground state for this system is a pure Fock 
state in configuration space, and therefore an eigenstate 
of the monopole operator and in fact of any operator of 
the form of Eq. jSJ, with very small corrections. This is 
of course a manifestation of the insulating nature of this 
state. 

The final figure is an overview picture. Figure [S] dis- 
plays the three moment-related frequencies in units of the 
lowest monopole frequency E2', the zeroth moment mg; 
and the quantum fluctuations in the local density at two 
points, all as functions of the coupling U and the trap fre- 
quency UJt ■ The data has been recorded along an oblique 
path through phase space along the line U = 20wt, where 
LUt has varied from the value 0.3 up to 4; this way, all the 
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FIG. 5: Properties of the system as functions of the coupling 
U and the trap strength ujt- The latter has been swept from 
LOt = 0.3 to ujt ~ 4, and and U has varied asU = 20a;f through 
the sweep. The topmost panel shows the three moment- 
related frequencies (from top to bottom luu), cj2o and uj2i) 
in units of the lowest monopole frequency £2- The middle 
panel shows the magnitude of the zeroth moment mo. The 
bottom panel plots the quantum fluctuations in local density, 
5n^ — (n^) — {n)^, at the point r = (solid line) and r — 2 
(dashed line) . The units are the dimensionless units employed 
throughout the paper. 

phases discussed above could be sampled. For small U 
and uJt, the system is superfluid and the density fluctu- 
ations are large. The moment-related frequencies stay 
close to each other and to the eigenfrequency i?2. The 
decrease of the density fluctuations in the centermost well 
signals the transition to a Mott-plateau state, and at the 
same time the frequencies begin to differ, but the zeroth 
moment is still above unity. The drop in mg comes at the 
same point as the drop in the fluctuations in the second 
off-center site; this marks the transition to the Mott in- 
sulating state. In this finite system, the Mott transition 



is gradual. 

Reference reports on an experiment that measures 
the response to excitation of trapped bosons in a three- 
dimensional optical lattice. It is found that for a shallow 
lattice, which corresponds to a relatively small U com- 
pared with the tunneling strength, one narrow frequency 
range is excited, just as expected for a trapped conden- 
sate. However, for larger U, a broad range of frequencies 
is excited with appreciable strength. This result, which 
is not easily interpreted in the framework of Bogoliubov 
theory, is in fact very similar to the present findings for 
a system with a Mott plateau. The drastic difference 
in size, dimensionality and method of excitation forbids 
quantitative comparison, but the qualitative similarities 
are there: we have seen in this paper that the rise and 
broadening of the frequencies is characteristic of being in 
the vicinity of the Mott-plateau regime. 

In conclusion, the strength function for the monopole 
and dipole excitations of the Hubbard model in a trap 
have been calculated, and with an outlook to treating 
larger systems, the moments of the strength function 
have been employed in order to extract information about 
the excitation strengths and frequencies. It is seen in 
both the moment calculations and the full strength func- 
tion how the finite system undergoes a transition be- 
tween superfluid, Mott plateau, and fully Mott insulating 
states. It is predicted that the experimental signature 
of a system containing a Mott plateau is the excitation 
of several frequencies with appreciable strength, while a 
strongly coupled system will not be excited at all because 
of its Mott insulating nature. Based on the present find- 
ings, it is suggested that the broad and strong resonance 
found in a recent experiment is due to a Mott insu- 
lating plateau in the system. 
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